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I. INTRODUCTION 

Recently ideas proposed by Lee and Wick [l|, |2| were used to extend the standard model so 
that it does not contain quadratic divergences in the Higgs mass Q|. Higher derivative kinetic 
terms are added for each of the standard model fields. They improve the convergence of 
Feynman diagrams and result in a theory where there are no quadratically divergent radiative 
corrections to the Higgs mass. The higher derivative terms give rise to propagators with new 
poles that are massive resonances. These Lee- Wick (LW) resonances have wrong-sign kinetic 
terms which naively give rise to unacceptable^nstabilities and violations of unitarity. Lee 
and Wick l|, |2| and Cutkowski et al. (CLOP) [4] proposed a way of defining the integrations 
that arise in Feynman diagrams so that the theory is unitary, Lorentz invariant, and free of 
instabilities. However, there is acausal behavior caused by the unusual location of poles in 
the propagators. Physically this acausality is associated with the future boundary condition 
needed to forbid the exponentially growing modes. As long as the masses and widths of 
the LW resonances are large enough, this acausality does not manifest itself on macroscopic 
scales and is not in conflict with scattering experiments. Various aspects of this model [sl,!^, 
.t. extension. M and of Lee-W.ck theone. gene^aliQAQQ have been explored 
in the recent literature. jCollider phenomenology |14J . 
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IgI . constraints from electroweak 



precision measurement [8|, [iTI, ll8|] and the cosmology of theories with higher derivatives pi 
have also been studied. 

In this paper we examine the high temperature behavior of Lee- Wick theories, including 
the LW standard model. In these theories the S-matrix can be calculated in perturbation 
theory using the prescriptions of Lee and Wick and CLOP. It is unclear whether a functional 
integral formulation of LW theory exists, so a computation of finite temperature effects 
solely based on the known S matrix is desired. The formalism of Dashen, Ma and Bernstein 
(DMB) [20] expresses the thermodynamic grand potential in terms of the S-matrix and we 
apply it to LW theories to deduce the pressure and energy density for these theories at 
finite temperature. Although previous analyses have argued that in scattering experiments 
no acausal effects persist to macroscopic scales, it is interesting to examine whether this is 
possible when multiple scattering effects play a role. This is the case for thermal equilibrium 
and we explore the propagation of sound waves in a gas consisting of ordinary and Lee- Wick 
particles. We find that in such a gas, at a large (but finite) temperature, sound waves 
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propagate at a speed less than light. 

In the next section we review scattering in a simple scalar Lee- Wick theory. Section 3 
uses the DMB formalism to calculate, in this toy model, the energy density and pressure 
at thermal equilibrium. At high temperatures T ^ M we find that the LW resonance 
contributes minus what an ordinary particle of mass M would. We use our results to 
conclude that as T — > oo the speed of sound approaches Cg — 1 from below. This gives 
in the hmit T — > oo a speed of sound equal to the speed of hght and is the largest value 
consistent with causal propagation of classical sound waves in the gas. Concluding remarks 
are made in Section 4. 



II. A TOY MODEL 



In this section we introduce a simple Lee- Wick theory with a single self-interacting real 
scalar field. In addition to the standard kinetic term there is a higher derivative term. The 
Lagrangian density is 

c = \d,^d^^ - ^,{d^r - \m^' - i#^ (1) 

so the propagator of in momentum space is given by 

For M ^ m, this propagator has poles at — and also at — M^. Thus, the 
propagator describes more than one degree of freedom. 

We can make these new degrees of freedom manifest in the Lagrangian density in a simple 
way. First, let us introduce an auxiliary scalar field cp, so that we can write the theory as 

C = ^djd''4> - \m^'' - 4>d'4> + Im'^' - ^g4>\ (3) 

Next, we define = 0-1-0. In terms of this variable, after integrating by parts, the Lagrangian 
density becomes 

£ = ld,<pd^<P - \d,~^d>^~^ + ^M202 _ lra\ct> - 0)2 - ^^(0 - 0)^ (4) 

In this form, it is clear that there are two kinds of scalar fields: a normal scalar field and 
a new field 0, which we will refer to as a LW field. The sign of the quadratic Lagrangian of 
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the LW field is opposite to the usual sign so one may worry about stability of the theory, 
even at the classical level. We will return to this point. If we ignore, for simplicity, the mass 
m, the propagator of is given by 

The LW field is associated with a non-positive definite norm on the Hilbert space, as indi- 
cated by the unusual sign of its propagator. Consequently, if this state were to be stable, 
unitarity of the S'-matrix would be violated. However, as emphasized by Lee and Wick, uni- 
tarity and Lorentz invariance can be preserved provided that (p may decay. This is natural 
in the theory described by Eq. (jlj) because is heavy and can decay into two particles. 

In the presence of the mass m, there is a mixing between the scalar field and the 
LW scalar 0. We can diagonalize this mixing without spoiling the diagonal form of the 
derivative terms by performing a hyperbolic rotation of the fields: = 0'cosh^ + 0'sinh^, 
= 0' sinh 6* + 0' cosh 6. This transformation diagonalizes the Lagrangian if 

, ^ -2mVM2 
tanh2^ = ttttt^- (6 

A solution for the angle 9 exists provided M > 2m. The Lagrangian density (jl]) describing 
the system becomes 

C = ^9,0'9>' - ^m'20'^ - ^9,0'9'^0' + Im"4>" - ^^'(0' - 0')', (7) 

where m' and M' are the masses of the diagonalized fields and g' = (cosh^ — sinh^)^5f. In 

what follows we assume that M ^ m, so that g' ^ g. 

Introducing the LW fields makes the physics of the theory clear. There are two fields; the 

heavy LW scalar decays to lighter scalars. At loop level, the presence of the heavier scalar 

improves the convergence of loop graphs at high energy consistent with our expectations 

from the higher derivative form of the theory. 

Loop corrections to the two point function of the LW field play a crucial role. Near 
= and at small g the (f) — (f) mixing can be neglected and the full LW propagator 

and its perturbative expansion are given by 

^"(P) = ^Zl^ + ^TTigi + ■ ■ • 

(8) 



p2 _ M2 - S(p2 
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Note that, unlike for ordinary scalars, there is a minus sign in front of the self-energy E(p^) in 
the denominator. The pole mass shift of the LW scalar coming from the radiative corrections 
is +S(M^). This sign is significant; for example, from a one-loop computation we see that 
the imaginary part of the self energy is 



ImS(/) = ^%^-4m^)./l-^. (9) 

Therefore the propagator develops a pole for Im(p^) > 0. In the narrow width approximation 
the propagator for the LW field is 

where 



This width differs in sign from widths of the usual unstable particles we encounter. Strictly 
speaking the propagator has an additional pole at p"^ = + iMT and a cut over the real 
axis; however, the effect of these two terms in the below calculation of the pressure and 
energy density at finite temperature cancel one another so we ignore them. 

Using the Lee and Wick and CLOP prescriptions the S-matrix in this theory is unitary 
and Lorentz invariant on the space of physical ordinary particles. The contour of integra- 
tion over p^ in this prescription does not lead to any instability, as one may naively guess 
from the negative sign in the width, but instead leads to apparently acausal behavior. This 
has been extensively discussed in the literature 

The LW resonance is unstable and therefore does not appear in the initial or final states 
of the S'-matrix. This is similar to the case of the VT-boson of the standard model, which does 
not appear in initial or final states of S-matrix elements because it is unstable^. Nonetheless, 
the resonance impacts the scattering of ordinary particles 0(pi) +0(^2) 4'{p'i) +0(^2); 
particularly near the kinematic point {pi +^2)^ = M'^. 

Writing S = 1 — iT, the T matrix can be computed using the standard Feynman tech- 
niques, modified appropriately for Lee- Wick theories. DMB introduce a closely related 



^ One difference is that there are no poles in amphtudes associated with the W^-boson two point function, 
just a cut that is represented as a pole in the narrow width approximation. However, for the Lee- Wick </> 
resonance there is actually a pair of negative residue poles in the complex plane in addition to the usual 
cut. 
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quantity T{E) that has two particle matrix elements 

(Pi, P2|T(i?)|pi, P2) = (27r)5(i? -E,- E^){2Tif5\^ - P)M{E), (12) 

where P = pi + p2, P = pi + P2 and M. is essentially the usual invariant matrix element. 
For center of mass energies near M the invariant matrix element is given (in the narrow 
resonance approximation) by^ 

Note that this differs from scattering via the exchange of an ordinary {i.e., not LW) resonance 
by an overall minus sign and the fact that F given by Eq. ( fTTl) is negative. 



III. THE PRESSURE AND ENERGY DENSITY IN THERMAL EQUILIBRIUM 

The grand partition function Vl at zero chemical potential is defined by 

n = ilnTre"^^, (14) 

P 

where (3 = 1/ kT and the trace is over all physical states in the theory. In our toy model the 
physical states are the particle states but not states that contain a LW resonance. From 
Q one can calculate the thermal equilibrium pressure p and energy density p in the usual 
fashion using formulas 

P = -y^ P=-^^ (15) 

where V is the volume of the system. 

DMB derive the following expression for the grand potential 

^ = ^0 - ^ / di^e-^-^ (rrASiEr^SiE)^ , (16) 

where c denotes that only connected diagrams are taken into account. In Eq. (fT6l) the 
S-matrix is given by S{E) = 1 — iT{E), Qq is the free particle grand potential and A is 
an operator that sums over permutations of the identical particles in the trace with the 
appropriate minus signs for fermions. Using the relation between S and T this becomes 



(3 J Am 



TrA I -.A [T{E)+r{E)^]+T{E)^^T{E) 



(17) 



^ There is an extra factor of 1/2 associated with identical particles that we have chosen to put in rather 
than in phase space integrations. 
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To evaluate Q in the toy model introduced in the previous section, we begin by evaluating 
the part of Q that comes from the contribution of two particle states. It is convenient to 
use the phase space relation 



(27r)32Ei7 (27r)32E2 J J {2'Kf2E[ J {2'nf2E'^ i ' ^ 
where the primed variables are the center of mass momenta and energies, uj = E[+ E2 and 



E = El + E2 = voJM-P^. In order to calculate the first term of the integral in Eq. (fT71) we 
notice that^ 

Tr^r(i?) = ^W{E), (19) 
so the expression we need to evaluate is 

Using (27r)^5'^(0) = V and the phase space relation f|T8l) gives 

(21) 

Now, recall that = (-E, P) is the total energy-momentum four-vector of the states in the 
trace while {u, 0) is the corresponding energy-momentum four- vector in the center of mass 
frame. They are related by a boost with a Lorentz gamma factor 7 = E/u. One can go 
from one set of variables to the other. Using the relation E"^ = u"^ + we have that 

^5iE - El - E2) = 5{uj -E[- E'2) (22) 
E 

and so the integrations over d'^Pid^p2 become the standard two body phase space integration. 
We arrive at the result 



Finally, we change from the variable to in all other places. Using 
and interchanging the order of the two integrations we get 



J dE ^ ' J (27r)3 J duo \ 87r\ 00^ ^ 



(25) 



^ For simplicity we drop the subscript c. 
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where 



1 



9' 



(26) 



2 cj2 - M2 + iMV ' 

The term with T{Ey gives the same contribution to Q but with Ai{uj) substituted by 
Ai*{uj). In calculating the second term of the integral in Eq. ( |T7I) we use the same methods 
as previously. We need to evaluate 



d 



TrT{Ey—T{E) = V j H 



n 



L (27r)32E, J f \ (27r)32E. 



{2ttY6\P - P) 



d' 



[{2tt)S{E - El - E2)M*{E)] — {27r)6{E - E^ - E2)M{E) 

oE L 



(27) 



where P = pi + P2, P = pi + P2 and the factor V{2tx)^5'^{P — P) came from the momentum 
delta functions in the definition (fT2|) of T{E). In the cm. frame we have 



T^nE)^-^T{E) = V 



d^P f d^P 



(27r)3 J (27r)3 

2 



(27r)^5-^(P 



d 



{2t:Y5{E~E,~E^)M*{E)\ — 



n ,/'fiA/ '(Pi + P2)^(2vr)^5(i? - - E,)M{E) 



V 



d^P 
(2^ 



1 / 4m2,^^^^^ 
Stt V oj^ 



\ (2vr)32K 
dE 



Stt V uj 



:28) 



Putting all this together and neglecting the mass m of the ordinary scalars, we arrive at the 
LW contribution to the grand potential in the form 



V f d^P 



LW 



P J (27r)3 



due-^-"^— 
Aixi 



d fMiuj) MHuj] 



duo \ Stt 



+ 



Stt 



+ 



MHuj) d Miuj] 



Stt duo Stt 



+. 



(29) 



where the ellipses denote the terms from summing over permutations which basically are 
multiple insertions of the two body state. Performing the differentiations and using the 
explicit formulas for M. and P this becomes 



V r d^p 



LW 



(3 J (27r)3 



duoe-^^^^ 



uoMT 



_TT (Cj2 - M2)2 + M2P2 



+ .... 



(30) 



Recall that for the LW resonance P is negative. The above formula is the same as one would 
get for scattering through an ordinary resonance except in that case P is positive. Therefore, 
in the narrow LW resonance approximation 

2 uoMT 



TT (a;2 - M2)2 + M2P2 



-S{uo-M). 



(31) 
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Hence, the contribution to the grand potential from the single LW resonance is 



single LW 



V 



(27r)3 



-/3VM2+P2 



(32) 



This is precisely what one would expect from a stable particle of mass M, except for the 
overall plus sign instead of a minus. Note that the narrow width approximation is valid 
provided the prefactor in the integral over u) is slowly varying. This will be the case provided 

(3\r\ < 1. 

Thus far we have only included two particle states in the calculation of the trace for the 
grand potential, with a resonant S matrix from Eq. (13). Contributions to the trace from 
states with more than two particle introduce multiple resonance amplitudes. This has two 
effects 221]: they modify the width, as expected for decays in a thermal bath, and they 
convert the exponential factor to the usual Bose-Einstein logarithm. In the narrow width 
approximation the modified width is still narrow, that is, proportional to g"^, and hence 
it still vanishes as g 0. The Bose-Einstein logarithm arises from considering multiple 
resonance graphs that are connected only because of the permutation operator A. This 
gives 

which is minus the contribution of a boson of mass M to the ideal gas grand potential. 
Note that this result is valid for arbitrarily large temperatures, which is not the case for 
non-elementary (composite) resonances 23|. 

Since in the narrow width approximation the LW resonance contributes minus what an 
ordinary scalar particle of mass M would, in our toy model the LW contribution to the 
energy density is 



Plw 



7i\kTy M\kTf 



30 



24 



+ 



while the contribution to the pressure is 
Plw 



+ 



(34) 



(35) 



90 24 127r 

Here the ellipses stand for terms of order ln(T) at most and are less important than those 
explicitly displayed when T ^ M. Adding these to the positive contributions from the 
ordinary scalar (whose mass m we neglect) 

TT^kT)^ TT^kry 



Pordinary 



30 



Pordinary 



90 



(36) 
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gives 



and 



P — Pordinary + PLW 



M^kTf 



24 



+ ... 



P — Pordinary + PLW 



M\kTy M^{kT) 



+ .... 



(37) 



(38) 



24 1 27r 

A similar analysis holds for theories with a left handed fermion and its LW partner. This 
is most easily seen in the higher derivative formulation of the theory. In the auxiliary field 
formulation one usually introduces left and right handed LW fermions but one of these is 
dependent on the other through the equations of motion. In that case^ 



and 



F 

Pun 



120 



24 



+ 



+ ..., 



(39) 



(40) 



360 24 

where, similarly to the boson case, the eUipses denote terms of order ln(T) at most. The 
absence of a term linear in T for the fermion pressure and density is as in the normal case. 
Adding the LW fermion contribution to the ordinary fermion energy density and pressure 



F ^ 77r\kTY 

Pordinary 



120 



p ^ 77r\kTY 

■^ordinary 



360 



(41) 



gives 



and 



F _ F F _ M^{kTf 

P — Pordinary + PlW " ^ ■ ■ ■ 



F_ F ^ F _M\kTf 

P — Pordinary + PlW " 



24 



+ .... 



(42) 



(43) 



From the above formulas for the pressure and energy density one can calculate the factor w 
in the equation of state p = to be 



w 



w 



1 

l + O 



2 M 

+ O 

TT kT 



ln(T) 



ln(T) 

J^2 



rp2 

for fermions 



for bosons. 



(44) 



In Figure 1 we plot w — p/ p for both cases as a function of kT/M. The value w — 1/3 
for small temperatures should not be puzzhng since for T — > the LW contribution is 



We take into account a factor of 2s + 1 in the formula for the grand partition function, where s is the spin. 
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FIG. 1: Factor w in the equation of state p = wp as a function of kT/M for fermions (solid) and 
bosons (dashed). 



suppressed by the Boltzmann weight factors and only the ordinary particles contribute to 
the grand potential. More interesting is the value w = 1 at high temperatures, which may 
have implications to the early universe cosmology. Cosmology with equation of state w = 1 
has been investigated in the context of holographic cosmology for a medium consisting of a 



black hole emulsion 



24 



25 



26 



27 



28|. 



The speed of sound Cg can be calculated using the formula 



dp dp 
df'df 



Taking into account higher order correction terms we arrive at 

'ln(T) 




(45) 



1 M 

V O 

2-K kT 



l + O 



ln(T) 
2^2 



2^2 

for fermions. 



for bosons, 



(46) 



Figure 2 shows the plot of Cs for fermions and bosons as a function of kT/M. The speed of 
sound increases from 1/VS at T = to 1 as T — oo, which is equal to the speed of light. 
Hence, the propagation of classical sound waves in the relativistic gas is causal. 
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FIG. 2: Speed of sound Cg = y^dp/dp as a function of kT/M for fermions (solid) and bosons 
(dashed). 

In the LW standard model there is a LW partner for every ordinary particle. Hence, at 
high temperatures above all the LW masses the pressure and energy density are approxi- 
mately proportional to M^T^, which implies the same equation of state w = 1. This gives 
a universe with density p ~ a(t)~^ and the scale factor a(t) oc t^^^. 

IV. CONCLUSIONS 

We have studied in this paper Lee- Wick theories at high temperature. Making use of 
the S"- matrix formulation of statistical mechanics presented in DMB [20i], we calculated 
the grand thermodynamical potential for a gas of Lee- Wick resonances in the boson and 
fermion case. We found that the contribution of Lee- Wick resonances to the energy density 
and pressure is negative for high temperatures. 

Next, we considered a gas of both ordinary and Lee- Wick particles. We found that for 
high temperatures (much greater than the mass of the resonance itself) the contributions 
of the Lee- Wick resonances to the pressure and energy density cancel against those of the 
normal particles at leading order in temperature {i.e., (9(T^)). We confirmed this for both 
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fermions and bosons. We found that the remaining 0{M'^T'^) contribution is positive and 
identical for the pressure and energy density. This led us to the equation of state w — 1 
for T ^ oo. In applications to big bang cosmology this yields a scale factor of the universe 

a{t) oc t^l'\ 

The quantity Cg = ^Aip/Vip corresponds to the speed of sound in the medium. We 
checked that is less than 1 in the whole temperature range, thus causality is not violated. 
It seems interesting to investigate the cosmological implications of the equation of state 
w — 1, especially for the propagation of fluctuations in the early universe. 
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